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Abstract. 

Most simulations of coronal mass ejections (CMEs) to date either focus on the interplanetary propagation of a giant plasma 
"blob" without paying too much attention to its origin and to the formation process or they focus on the complex evolution of 
the coronal magnetic field due to (sub-)photospheric motions which result in an eruption. Here, we present global simulations 
of CMEs where coronal motions are used to produce a realistic evolution of the coronal magnetic field and cause an eruption. 
We focus on active region 10069, which produced a number of eruptions in late August 2002, including the August 24, 2002 
CME - a fast (~ 2000 km s _1 ) eruption originating from W81-, as well as a slower eruption on August 22, 2002 (originating 
from W62). Using a three-dimensional magneto-hydrodynamic (MHD) simulation of these ejections with the Space Weather 
Modeling Framework (SWMF), we show how a realistic initiation mechanism enables us to study the deflection of the CME 
in the corona and in the helio sphere. Reconnection of the erupting magnetic field with that of neighboring streamers and active 
regions modify the solar connectivity of the field lines connecting to Earth and change the expected solar energetic particle 
fluxes. Comparing the results at 1 AU of our simulations with in situ observations by the ACE spacecraft, we propose an 
alternate solar origin for the shock wave observed at LI on August 26. 
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1. INTRODUCTION 

1.1. Active Region 10069 and August 
Eruptions 

On August 24, 2002, active region (AR) 10069 was 
near the western limb of the Sun (W81) when it was the 
source of a fast (~ 1,900 km s _1 ) and wide coronal mass 
ejection (CME). This event has been well studied due to 
extensive observations remotely by LASCO and UVCS 
[1] as well in-situ by the Wind and ACE spacecraft. Most 
importantly, it was associated with an intense Solar En- 
ergetic Particle (SEP) event [e.g., see 2, 3]. The eruption 
was preceded 2 days earlier by a similar, but twice slower 
CME, which was also associated with a large SEP event 
but is not thought to have been associated with a shock 
wave at 1 AU. Particularly noteworthy is the fact that, at 
the start of the August 24 event, the proton flux at ener- 
gies above 10 MeV was still ten times above the intensity 
before the August 22 CME. 

Western limb events such as the August 24, 2002 CME 
present a number of challenges for space weather predic- 
tion. Due to the Parker spiral, the Earth is magnetically 
connected with regions on the solar surface at around 
W55 (for a 400 km s _1 wind). Therefore, SEP events 
are preferentially associated with western-limb events. 
Events such as the August 24, 2002 CME present ad- 



ditional challenges since the SEP arrival time at Earth 
corresponds to a particle release height of less than 5 R Q 
[2]. To explain the observations, a shock wave must be 
formed low in the corona at the flanks of the CME or the 
magnetic field line connecting the Earth to the solar sur- 
face must significantly diverge from the nominal Parker 
spiral. Such large differences of up to 30° between flare 
sites and the footprint tracing of the magnetic field lines 
connecting the Earth to the Sun have been reported be- 
fore [4]. 

Western limb events can be associated with ejecta and 
more often with shock waves at 1 AU [5]. It is believed 
that the August 24, 2002 CME was associated with a 
shock wave at 1 AU on August 26. This fact, again, 
seems to imply either a very large span of the shock 
wave, a large deflection of the CME, or a combination 
of both. However, in Lugaz et al. [6], we studied numeri- 
cally the August 24, 2002 eruption, and our model could 
not explain the expansion and/or deflection required to 
associate the shock at 1 AU to the eruption on August 
24. In this article, we investigate i) why the August 22 
CME, which originated from the same region when it 
was closer to the disk center, was not associated with a 
shock at 1 AU, and ii) the importance of this preceding 
ejection on the transport of SEPs accelerated by the Au- 
gust 24 CME. 
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FIGURE 1. Top: Magnetic topology of the solar corona in 
the vicinity of AR 10069 on August 22-24, 2002. Bottom: 
Snapshot of the eruption at the end of the shearing phase (t = 1 
hour) showing magnetic field lines color-coded with the radial 
velocity. The reconnection between the erupting flux from the 
dipole and the open magnetic field line happens at eNP and is 
highlighted with the white circle. 

1.2. Solar Wind and Solar Eruption Models 

Our simulation setup is similar to that discussed in 
more details in Roussev et al. [7], Jacobs et al. [8] and 
Lugaz et al. [6]. In our numerical model, the steady-state 
solar corona and solar wind are constructed following 
the methodology of Roussev et al. [9]. The initial condi- 
tion for the coronal magnetic field is calculated by means 
of potential field extrapolation with boundary conditions 
for the radial magnetic field at the Sun provided by full- 
disk SoHO/MDI observations. The plasma parameters 
are prescribed in an ad-hoc manner, through a variable 
polytropic index, in order to mimic the physical proper- 
ties of streamers and coronal holes once a (non-potential) 
steady state is reached. 



To the initial magnetic field constrained by MDI data, 
we superimpose newly emerged magnetic flux given by 
the dipolar magnetic field of two point charges, whose 
magnitude is chosen so that the peak value of the radial 
magnetic field at the solar surface is about 47 Gauss. 
To initiate the eruption, we use a similar strategy as 
described in Roussev et al. [7]. To summarize, once 
the steady-state is reached at t = 0, the two subsurface 
magnetic charges are moved apart quasi- steadily up to 
t = ts = 60 min with a speed which is ramped up in t$/3 
to V q = 33 km s _1 ; the charge motion is stopped at t = ts. 
The two charges are moved approximatively parallel to 
the polarity inversion line. In order to provide sufficient 
buildup of magnetic energy, accompanying horizontal 
boundary motions are imposed as well. 

2. RESULTS 
2.1. Loss of Equilibrium 

One of the main results of the work in Roussev et al. 
[7] was to recognize the importance of the pre-existing 
magnetic topology in the initiation of the eruption. This 
is first and foremost because reconnection at the pre- 
existing null points (NP) and quasi- separatrix layers 
(QSLs) enables the sheared and energized magnetic flux 
of the dipole to erupt. There are four main pre-existing 
topological features important to understand these two 
eruptions (see top panel of Figure 1): a northern NP 
(nNP) between ARs 10067 and 10069 and a QSL be- 
tween ARs 10066, 10068 and 10069 for the August 24, 
2002 eruption and two NPs for the August 22 eruption a 
southern one (sNP) between ARs 10069 and 10079 and 
an eastern one (eNP) between ARs 10079 and 10083. 

For the August 24 eruption, more details can be found 
in Lugaz et al. [6]. As a summary, the eruption is trig- 
gered when expanding field lines reconnect with the 
overlaying field through the QSL. As some of these field 
lines expand further, they reconnect though the nNP and 
some of them open up. The main motion of the erupting 
flux is radially above the initial position of the QSL but 
reconnection through the nNP enables the expansion of 
the CME towards solar east, i.e. toward the Earth's direc- 
tion. 

For the August 22 eruption, the EUV images show 
that the flare and CME occurred around the eastern side 
of AR 10069. Consequently, the two magnetic charges 
are put on the eastern side of AR 10069, approxima- 
tively midway between nNP and sNP. As the charges are 
moved, the initial reconnection, triggering the eruption, 
occurs at sNP and the eruption moves radially above this 
NP. Therefore, the initial direction of the August 22 erup- 
tion is to the south-east of the direction of the August 24 



eruption, as also appears on the coronagraphic images. 
As the erupting field lines expand, they start to reconnect 
with open field lines at the eNP, resulting in magnetic 
connection between the source region of the eruption and 
the Earth (see bottom panel of Figure 1). 



2.2. Magnetic Connection to Earth 

At the start of the August 22 eruption, Earth was 
magnetically connected to the Sun through the eNP. The 
structure of this NP is such that the field lines in a 
cone 10° around Earth have footprints spanning from 
N20W50 to SI 5E 10 on August 22. This explains why 
this S07W62 eruptions was very well connected to Earth 
and could result in a large SEP event. The flare site was 
only 20° away from open field lines connected to the 
vicinity of Earth. During the second half of the shearing 
phase, the erupting flux reconnects with some of these 
open magnetic field lines through the eNP (see bottom 
panel of Figure 1). This creates an opportunity for flare 
energetic particles to be further accelerated through the 
shock wave, creating a mix of flare- and shock-related 
particles (or more precisely of shocked superthermals 
and shocked flare particles). 

When we consider the magnetic connectivity between 
the source region of the ejection and the Earth at the start 
of the second eruption, we cannot ignore the presence 
of the previous ejection, which is at a distance of almost 
0.95 AU (see bottom panel of Figure 2). Field lines con- 
necting Earth and the Sun are still passing through the 
eNP, but they are modified due to the presence of a pres- 
sure disturbance associated with the previous eruption. 
We find that the length of the field line passing through 
Earth is significantly shorter when taking into account 
the August 22 CME compared to when it is calculated 
from the steady-sate solution (1.06 AU compared to 1.38 
AU). This is because the field line goes through faster 
wind associated with the CME without being "hindered" 
by the closed field lines of the ejecta itself. This shorter 
length should be taken into account when deriving the 
release time of the particle as well as the height at which 
the shock wave associated with the August 24 eruption 
has formed. It should be noted that our steady-state solar 
wind is too fast compared to observations at 1 AU, and 
the precise values given here are therefore not reliable. 



2.3. Heliospheric Evolution 

A shock wave was detected at Earth by the ACE space- 
craft about 58 hours after the August 24 eruption. No 
shock wave was detected during the period August 19- 
25. One of the goals of our investigation is to deter- 
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FIGURE 2. Top: Erupting magnetic field lines at the end of 
the shearing phase (t = 1 hour) color-coded with velocity. Note 
the absence of twist in the magnetic field. Bottom: Aug. 22 
CME at the start time of the August 24 CME. The view is 
in the ecliptic plane, the large black disk represents the Sun, 
and the small one near (-140,-170) represents the Earth. The 
black circle has a radius of 1 AU. The contours show the density 
scaled by 1/r 2 . 



mine whether it is possible that the August 22 erup- 
tion (W62) does not result in a shock wave detectable at 
Earth, whereas the August 24 eruption (W81) is associ- 
ated with a shock detectable at Earth. As shown in Lugaz 
et al. [6], in our simulation, the shock wave associated 
with the August 24 eruption does not extend azimuthally 
all the way to Earth, even though it has an angular span 
of about 120°. We also want to investigate whether or 
not the August 22 CME could be the cause of this shock 
wave. If it is not, the presence of a disturbed solar wind 
ahead of the August 24 CME, may explain at the same 
time the 58 hours transit time and the strong deceleration 
of the shock driven by the August 24 CME. As measured 
at 1 AU, the downstream region of the shock wave had 



a speed of about 450 km s _1 with a jump in speed of ~ 
50 km s -1 [2]. Looking at the August 22 CME and its as- 
sociated shock wave in our simulation (bottom panel of 
Figure 2), we find the opposite problem: the shock is de- 
tectable at 1 AU, 50 hours after the ejection, which shows 
that our CME model has too much azimuthal expansion. 
We should note that the solar wind speed is about 150 
km s _1 too large before the arrival of the shock wave. 
Although this might partially explain the early arrival 
time, it is not enough to double it to match the obser- 
vations. It is also worth mentioning that the Empirical 
Shock Arrival (ESA) model [10, 11] for a CME whose 
initial speed is 1000 km s _1 (resp. 1900 km s _1 ) gives 
a transit about 52 hours (resp. 24 hours). One of the few 
ways to associate the shock wave observed at Earth with 
one of these two CMEs, is to explain it as the flank of the 
shock wave driven by the August 24 CME, if its initial 
speed is 900 km s _1 . This would imply a very wide and 
very curved shock front with a nose initially traveling at 
2000 km s" 1 and the flanks 90° away at 900 km s" 1 . In 
general, our simulations fail totally in explaining the ar- 
rival of a shock wave at Earth on April 26 associated with 
one of these two eruptions. The August 24 CME is not 
deflected enough and the August 22 CME cannot decel- 
erate enough to arrive at Earth in over 4 days. Therefore, 
we look back at the CME list on August 22-24 to find a 
possible candidate. We believe the disk-center CME/flare 
(S4E22) on August 23 at 0548 UT is a good candidate for 
the driver of the shock observed on August 26. It appears 
to be associated with a partial halo whose speed is about 
700-1000 km s- 1 . A 700 km s" 1 CME is expected (ESA) 
to drive a shock which would arrive at Earth in about 74 
hours, in good agreement with what is observed. 

3. DISCUSSIONS AND CONCLUSION 

We have performed a Sun-to-Earth simulation of the Au- 
gust 22, 2002 CME event with a realistic CME initia- 
tion mechanism [7], following a similar simulation of 
the faster and more studied August 24, 2002 CME [6]. 
We have found that, although the August 22 CME was 
not detected at Earth, it is likely that the presence of its 
associated shock wave and ejecta in the heliosphere sig- 
nificantly modified the length of the field lines connected 
to Earth at the start of the August 24 CME. This, in turn, 
must be taken into account when deriving the SEP re- 
lease time and the shock formation height. We have also 
investigated the origin of the shock wave observed in situ 
at 1 AU on August 26 near 14UT. Our simulations are 
not able to reproduce the correct transit time and hit/miss 
from these two eruptions: the August 24, 2002 ejection 
is not wide enough to hit Earth in our simulation and 
the August 22, 2002 CME is predicted to hit Earth but 
1.5 days earlier than what was detected. It is possible 



that complex interaction in the heliosphere between these 
two eruptions and other eruptions in the period from Au- 
gust 22-24, 2002 modified the transit time of one of the 
shocks to explain the observations at Earth [as described 
in 12]. A disk-centered eruption on August 23 may play 
an important role in the formation of the detected shock. 

This study has shown the advantages and the limita- 
tions of modeling to study complex events such as fast, 
wide CMEs from the limb. The simulation can give us 
an estimate of the length of the field line connected to 
the Earth which can be, for example, compared to the 
ones computed from anisotropic measurements [13] and 
must be used to derive SEP release time and the height 
at which the CME-driven shock wave was formed. Here, 
we have found that the length of the field line decreases 
by about 20%. On the other hand, the model is not able to 
help to explain the arrival of a shock wave 2.5 days after 
the August 24 CME at a location about 80° east of the 
flare location. Although such Sun-Earth events originat- 
ing from the western limb are rare, it will be important 
to better understand them in the future in order to have 
more reliable space weather prediction. 
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